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Abstract 

We point out how, in certain models of New Physics, the same combination of couplings occurs 
in the amplitudes for both D^-D^ mixing and the rare decays l+l' . If the New Physics 

dominates and is responsible for the observed mixing, then a very simple correlation exists between 
the magnitudes of each; in fact the rates for the decay £~^i~ are completely fixed by the 

mixing. Observation of i^£~ in excess of the Standard Model prediction could identify New 

Physics contributions to D^-D^ mixing. 
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I. INTRODUCTION 



Following many years of effort, there is now indisputable experimental evidence for D^-D^ 
mixing. The current values (the HFAG 'no CPV-allowed' fit of the mixing parameters 
are 

- ^ = 0.0100l°:r and Vb ^ ^ = 0.0076t',Zl ■ (1) 

J- D D 

These show that (i) charm mixing occurs at about the percent level, (ii) xd, ?/d are compara- 
ble in magnitude and (iii) the signs of xd and ?/d are positive (although a direct measurement 
of the sign of xy) is yet to be made). 

While it is quite likely that the observed mixing amplitude is dominated by the Standard 
Model contributions, the exact predictions are quite difficult.^ There are several reasons 



for this 



aaa. 



For example, in the "short distance" approach [5i| at leading order in the 
Operator Product Expansion (OPE) formalism (operators of dimension D = 6), the individ- 
ual diagrams are CKM-suppressed to the level 0{X^) (A ~ 0.22 is the familiar Wolfenstein 
parameter), hinting that the observed charm mixing is a simple consequence of CKM struc- 
ture. This is, however, not correct because severe cancellations between diagrams (even 
through 0{as)) greatly reduce the D = Q mixing to 0(10"^) 0, 4]. As for higher [D > 6) 
orders in OPE, it is true that certain enhanced contributions have been identified 0, Q] , but 
a definitive evaluation is lacking due to the large number of D > 6 operators and the inabil- 
ity to determine their matrix elements. A promising alternative approach which involves a 
hadron-level description may be able to account for the observed magnitude of |/d and 
Xd, but predicts their relative sign to be opposite. It is fair to say that this is probably not 
the final word on the SM analysis. 

n 

Given the uncertain status of the SM description, it would be tempting but premature [9|] 
to attribute the observed xb to New Physics.^ But clearly, the possibility that NP makes 
a significant or even dominant contribution to the observed mixing is open. A recent com- 
prehensive treatment of NP models shows that a large number of such models can 



^ Henceforth, we will make frequent use of the abbreviations SM for Standard Model and NP for New 
Physics. 

^ We will focus on xb in this paper. Not only does the SM estimate for yu work reasonably well when long 
distance effects are included [10|, but it has also been shown that NP effects are too small to have any 
significant impact [HI. [iR [ist . 
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accommodate a value of xd at the per cent level. This encourages us to further explore the 
NP option. In particular, New Physics could affect charm-related processes beyond mixing, 
such as rare decays |l5]. Of special interest are the — > decays. At present, there 



are only the upper limits 



16 



3,Q, 



19| 



Bdo^^+^,- < 1.3 X 10"^ Boo^^+e- < 1.2 X 10-^ and 

BDO^f,±e^ < 8.1 X 10-^ , (2) 

all at CL=90%. Such branching fractions place bounds on possible NP couplings, which can 
be compared with that obtained from D^-D^ mixing. In this paper we study the impact of 
NP on the combined system of D^-D^ mixing and the rare decays It should 

be stressed that the SM rate for the decay mode, Boo^^+^- ^ 3 x 10^^^ can be estimated 
fairly reliably even upon accounting for the effect of long distance enhancement [2^ . This 
smallness of the SM signal makes it easier for NP contributions to stand out. In this paper, 
we point out how, in certain NP models, the same couplings occur in the amplitudes for 
both D^-D^ mixing and —* i^£^ decay. If the NP effects are significant in mixing, then 
a correlation will exist between the magnitudes of each. In fact the correlation can be very 
simple and striking, with the branching fraction B£)0^i+£- being proportional to the mixing 
parameter xb- 

For each NP model considered in this paper, we shall make the simplifying assumption 
that the NP dominates D^-D^ mixing and then derive the correlated branching fraction 
that is then predicted. Obviously, if NP does not dominate, all our results for the branching 
fractions become upper bounds. We shall stress the general issue of which conditions allow 
for such correlations and give specific examples. Finally, even if the number of parameters 
in a given NP model is too large to give a unique prediction for Boo^e+i- in terms of xd 
{e.g. Z' models, etc), we show in Sects. Ill B-D how it is possible to estimate the scale of 
Bdo^i+i- by using the value of xd as input. 



II. EFFECTIVE LAGRANGIANS 



Heavy particles present in NP models are not produced in final states of charm quark 
decays. Yet, effects generated by exchanges of these new particles can be accounted for 
in effective operators built out of the SM degrees of freedom. That is, by integrating out 
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degrees of freedom associated with new interactions at a heavy scale M, we obtain an 
effective hamiltonian written in the form of a series of operators of increasing dimension. 
Here, we restrict our attention to the leading order operators, of dimension D = 6. For both 
jjo _ j^o ]2iixing and —>■ decays, the complete basis of effective operators is known 
and is expressed most conveniently in terms of chiral quark fields, 

{f\nNP\t) = Gj2Q{fi) (/|g,K)(/i) , (3) 

i=l 

where the prefactor G has the dimension of inverse-squared mass, the Cj are dimensionless 
Wilson coefficients, and the Qi are the effective operators of dimension six. Throughout, our 
convention for defining chiral projections for a field q{x) will be qL,R{,x) = (1 ± 7 .0q (y)/2. 



For AC = 2 processes, there are eight effective operators that can contribute 14j. l2l|. 



(4) 

= {ulCr) (urCl) , Qj = {ulCr) {ulCr) , 

Q4 = (urCl) (urCl) , Qs = iuL(T^„CR) (uLfT^^CR) . 

These operators are generated at the scale M where the NP is integrated out. A non-trivial 
operator mixing then occurs via renormalization group running of these operators between 
the heavy scale M and the light scale /i at which hadronic matrix elements are computed. 

All possible NP contributions to c — > ui~^i~ can be similarly summarized. In this case, 
however, there are now ten operators, 

Ql = (JlI/l) (uLrCL) , Q4 = (44) (urCl) , 

Q2 = (47/l) (urYcr) , Q, = (iRcr^jL) {urct^'cl) , (5) 
Q3 = (urCl) , 

with five additional operators Qq, . . . , Qio being obtained respectively from those in Eq. ([5]) 
by the substitutions L ^ R and R ^ L. The corresponding Wilson coefficients will be 
denoted as Ci{fi). It is worth noting that only eight operators contribute to , as 

{i+i-\Q,\D^) = {i+i-\Qio\D'') = 0. 

To obtain a general expression for xd as implied by the effective Hamiltonian of Eq. Q, 
we evaluate the D^-to-D^ matrix element in the modified vacuum saturation approximation 



of Appendix A and work at the light scale /i = rric 



xb = G 
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(6) 



where we have taken A^c = 3, we remind the reader that rj is discussed in Appendix A and 
the prefactor G defines the scale at which NP is integrated out. To use this expression one 
must relate the light-scale coefficients {Ci(mc)} to their heavy-scale counterparts {Cj(M)} 
in terms of the RG-running factors given in Appendix A. 

The rare decays i'^l" and D'^ fi~^e~ are treated analogously. To the decay 

amplitude 

M = u{p_,s.)[A + B-f5]v{p+,s+) , (7) 
are associated the branching fractions 



Bdo- 



Bdo- 



D 
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Ml 
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\A\^ + \B\^ 



(8) 



where electron mass has been neglected in the latter expression. Any NP contribution 
described by the operators of Eq. gives for the amplitudes A and B, 



\A\ = G 



\B\ = G 



Arrir 



C3-8 + C4-9 



D 



Ml 
rrir. 



{Ca-3 + Cg-s 



(9) 



with Cj_fc = d — Gk- In general, one cannot predict the rare decay rate by knowing just 
the mixing rate, even if both xp and Bpa^^+i- are dominated by a given NP contribution. 
We shall see, however, that this is possible for a restricted subset of NP models. 



III. NP MODELS WITH TREE-LEVEL AMPLITUDES 

This is the most obvious situation for producing a correlation between mixing and decay 
because there is a factorization between the initial and final interaction vertices. In the 
following, it will be convenient to consider separately the propagation of a spin-1 boson V 
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and of a spin-0 boson S as the intermediate particle in the tree-level amplitudes. The bosons 
V and S can be of either parity. 

Spin-1 Boson V: Assuming that the spin-1 particle V has flavor-changing couplings and 
keeping all the operators in the effective Lagrangian up to dimension 5, the most general 
Lagrangian can be written as 

nv = nl^^^ + n^ , (lo) 

where the quark part Ti.y^^'^ is 

n^^^^ = gviULl^^CLV^ + gv2URi^cnV^' + gvzUL(y^..CRV^'' + gv^URa^.CLV^'' (11) 

and the part that describes interactions of V with leptons Tiy is 

= gW^Ll^^^iV^ + gW^Ri/RV^' + g'y~^L(y,u^RV^'' + gW^RCJ^jLV^r (12) 

Here is the vector field and V^y = d^V^ — d^V^ + ... is the field-strength tensor for V^. For 
this study it is not important whether the field V corresponds to an abelian or non-abelian 
gauge symmetry group. 

In order to see the leading contribution to D mixing from Eq. (|TT1) . let us consider a 
correlator, 

S^(g2) = (_^) j e*(5-^')-"(D°(p)|T {n^^=\x)n^^=\0)} |/^°(p)). (13) 
This correlator is related to the mass and lifetime differences of a D-meson as 

Ez)(M2) = 2Md (aMd - ^ATd) . (14) 



Inserting Eq. fllip for Ti,^'^ ^, we obtain 

SD(g') = H) P I d^x e'^''-P>''{0\T{V'^{x)ViO)} |0)(S"(p)| g^, uliM^) UL^.CLiO) 

+ gVigV2 Ul1^Cl{x) Ur'J„Cr{0) +gvigV2 Ur'J^Cr{x) ULluCLiO) 

+ gv2 URlf.CR{x) urJuCr{0) + 0{l/Mv) \D\p)) , (15) 

where 0{1/My) denotes terms additionally suppressed by powers of 1/My and (3 ~ 0{1) 
denotes all relevant traces over group indices associated with V^, with (3 = 1 for an Abelian 
symmetry group. The leading-order D = Q contribution is found by expanding the vec- 
tor boson propagator in the large My-limit and then performing the resulting elementary 



integral, 



/5 



D) 



Ml 



{D\p)\ 9'v,Qi + 2gvi9v2Q2 + 9t^2Q6 • 



(16) 



Taking into account RG-running between the heavy scale My and the light scale fi = rric at 
which the matrix elements are computed, we obtain a subcase of the general Eq. (I6l), 



X 



(V) 
D 



2MITd 



(Ci(m,) +C6(mJ) - 



■1 




C2(m,) + 


■ 1 






.2 






.12 







(17) 



where the superscript on xj^'' denotes propagation of a vector boson in the tree amplitude. 
The Wilson coefficients evaluated at scale /i = are 



Ci (rric) = r{mc,My) g^^ 



9vi9v2 



C2(mc) = 2 r(mc, Mvy^'^gvigv2 , C6(mc) = r(mc. My) 5-^2 



(18) 



Similar calculations can be performed for the i'^i decay. The effective Hamilto- 

nian in this case is 



M2 



gvig'viQi + gvig'v2Q7 + g'v\gv2Q2 + gv2gv2Qt 



(19) 



which leads to the branching fraction. 



^(V) 



1 - ^TTT^ {gvi - gv2f {g'vi - g'v2) 



(20) 



" 327rM^rfl ^ Ml 

Clearly, Eqs. f|T7|) and fl20|) can be related to each other only for a specific set of NP models. 
We shall consider those shortly. 

Spin-0 Boson S: Analogous procedures can be followed if now the FCNC is generated 
by quarks interacting with spin-0 particles. Again, assuming that the spin-0 particle S has 
fiavor-changing couplings and keeping all the operators in the effective Hamiltonian up to 
dimension five, we can write the most general Hamiltonian as 



lis = rLg + rts , 
where the quark FCNC part is given by 

T^^cNc ^ gsiULCRS + gs2URCLS + gssULlfMCid'^S + gsiUR'-f^CRd^'S 



(21) 



(22) 



and the part that is responsible for the interactions of 5* with leptons is 

= 9sMrS + 9's2^dLS + g's^^Ll^iLd^S + g'sJul^lRd^S . (23) 

Inserting this Hamiltonian into the correlator Eq. (1131) and performing steps similar to the 
spin-one case leads to 



Si,(M: 



1 



Ml 



(DO(p)I gl.Q-r + 2gsigs2Q3 + g^Q^ \D\p) . 



(24) 



Evaluation at scale fi = rric gives 



X 



(S) 
D 



2TdMI 



1 ri 

h - 

12 2, 
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Csim^) - Y^(C4(me) + Crim^)) + 77(^5 (mj + Csim,)) 



(25) 



with the Wilson coefficients defined as 



Clinic 



-2r(mc, Ms) ^ gsigs2 

r+{m„Ms)+ [- + 



\2 V24T 
1 



8^241 

'1 



2 V24T. 

r+(mc, Ms) - r_(mc, Ms)] gl2 

o \ (I ° ' 

r+{m,,Ms)+ - + 



r_{mc,Ms) 



9I2 



(26) 



2 ^/24T 



2 V24T 



r_(mc,M5) 



^51 



1 



[r+{m„Ms)-r_{m„Ms)]gl^ 



8V24T 

where for notational simplicity we have defined r± = 7'(i±v^24T)/6 (^(^ygq^ (Al)). 
The effective Hamiltonian for the decay is 

1 



n 



(S) 



Ml 



gsig'siQ^ + gsig's2Q8 + g'si9s2Qz + gs2g's2Q^ 



(27) 



and from this, it follows that the branching fraction is 



B 



,(S) 



D 



1287rm2M|rB \ M 



1 - 



4m« 



(5'5i - gs2)' 



ia'si + %2)^ 1 - 



Ml 



i2i 



Note that if the spin-0 particle 5* only has scalar FCNC couplings, i.e. gsi = gs2, no 
contribution to — > i^i^ branching ratio is generated at the tree level; the non-zero 
contribution to rare decays is produced at one-loop level. This follows from the pseudoscalar 
nature of the Z)-meson. 
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If there are not one but several particles mediating those processes (assuming that they 
all couple to quarks and leptons), the above generic Lagrangians would need to be modified. 
For example, in the spin-0 case, one would have to replace gi S with J2k 9ikSkF{k), where 
F{k) is a numerical factor, as Sk are the mediating fields. For example, in models with extra 
dimensions the factor F{k) would be related to Kaluza-Klein decompositions of bosons living 
in the bulk. Similar corrections have to be performed in the case of a bulk spin-1 boson. 

Below, we consider generic models where the correlations between the — mixing 
rates and — > rare decays can be found. 



A. Heavy Vector-like Quarks: Q = +2/3 Singlet Quark 



Scenarios with heavy quarks beyond the three generations are severely constrained ex- 
perimentally, if those quarks have chiral couplings. We thus examine the case where the 
heavy quarks are SU{2)l singlets (so-called vector-like quarks) [2^. Here, we consider the 
charge assignment Q = +2/3 for the heavy quark and then Q = —1/3 in the next Section. 



23 



2J] in which the Stan- 



Weak isosinglets with Q = +2/3 occur in Little Higgs theories 
dard Model Higgs boson is a pseudo-Goldstone boson, and the heavy iso-singlet T quark 
cancels the quadratic divergences generated by the top quark when performing quantum 
corrections to the mass of the Higgs boson. Weak isosinglets with Q = —1/3 appear in Eq 



GUTs 
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261], with one for each of the three generations {D, S, and B). 



The presence of such quarks violates the Glashow-Weinberg-Paschos naturalness condi- 
tions for neutral currents [27]. Since their electroweak quantum number assignments are 
different than those for the SM fermions, flavor changing neutral current interactions are 
generated in the left-handed up-quark sector. Thus, in addition to the charged current 
interaction 



(ch) 
int 



9 

V2 



there are also FCNC couplings with the boson 22(] 



(ntl) 



9 



int 



(29) 



(30) 



2v^cos^^ 

Here, g is the SM SU{2) gauge coupling and V^i is a 4 x 3 mixing matrix with a running 
over 1 ^ 4, z = 1 ^ 3, with the CKM matrix comprising the flrst 3x3 block. 
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FIG. 1: (a) D^-D^ Mixing, (b) fi+ij-. 

D^-D^ Mixing: In this case, a tree-level contribution to AM^i is generated from Z^-exchange 
as shown in Fig. [H This is represented by an effective hamiltonian at the scale Mz as 

= Scos^e.Mf^'^ = ^f^^ ' ^^^^ 

where from unitarity, 

Kc^-iv:,v,,+v:,v,s + v:,v,,) . (32) 

Thus, we find 

(+2/3) _ 2GFXlJl,MDBDr{m,,Mz) 

Using r{mc,Mz) = 0.778 and demanding that the NP contribution is responsible for the 
observed mixing value yields A^c = 2.39 x 10^"^. 

fj^^fi~ Decay: For this case, the leptons have SM couplings to the Z'^. We then have 

Aoo^e+e- = B^o^i+i- = Xuc — - ■ (34) 

Restricting our attention to the fi^fi~ final state because the decay amplitude is proportional 
to lepton mass, we find 



A^c ^ 8itBdo^u+, ° 



2 



1 - 



Mr 



2 1-1/2 



(35) 



From the branching fraction bound of Eq. ([2]), we obtain A^c < 4.17 x 10"^, which is much 
less restrictive than the value from D° mixing. 

Combining the Mixing and Decay Relations: A correlation will exist in this case because 
the coupling between Z° and the lepton pair is known from the Standard Model. Thus, 
if we assume that all the D meson mixing comes from the Q = +2/3 heavy quark {i.e. 
_^(+2/3) _ ^^-j^ then we can remove the dependence on the NP parameter A„c and predict 
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Boo^fi+fi- in terms of xd, 



3V2 

64^ Bur{m,, Mz) 



4.3 X 10"% < 4.3 X 10^ 



4mf, 
Mn 



1/2 



-11 



(36) 



B. New Gauge Boson Z' 



New heavy neutral gauge bosons can exist in a variety of NP models [28|] . In these scenar- 
ios, there are in general five parameters that describe the two processes under consideration 
here, namely gz'i, g'z'iy 9'z'2^ Mz', where the coupling constants are defined as 
in Eqs. ( ITTl [T2l) by substiuting V Z'. There are, of course, many ways to reduce this 
number. In the following, let us assume that Z' couples only to left-handed quarks and has 
SM-like diagonal couplings to leptons, 

9z'2 = 0, gz'i = ^ -- + sm Ow , gz'2 = n — ' i^V 

cos Uw \ 2 ) COSt^VK 

where g is again the SM SU{2) gauge coupling. This procedure reduces the number of 
unknowns to two, gz'i and Mz'- Note that for purely vector couplings of a Z' to leptons, 
i.e. g'zii = g'z'2 no contributions are generated for —>■ ^~ due to conservation of vector 
current. 

D^-D^ Mixing: The contribution of the Z' model to mixing is given by Eq. (|T71) . 

(z') _ pDMDBDr{mc,Mz') gh, 

" 3f^ Ml, • ^^^^ 

For the very slowly varying RG factor, we have taken r{mc, Mz') = 0.71, which is typical 
of values for a Z' mass in the TeV range. From Eq. (l38l) . we obtain the bound Mz'/gz'i > 
1.7 X 10*^ GeV. 

fi~^fi~ Decay: In this model, the contribution to the rare decay branching fraction can 
be written in the form, 



r(z') _ 



GpflmlMD 



Ami all Ml ^gg^ 



Besides the gl^^/Ml dependence which appears in the above D mixing relation of Eq. (138|) . 
there is now an additional factor of M|/M|/. The bound obtained from Eq. ([2]) implies 
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the restriction Mz' j gi\ > 8.7 10^ GeV, which is weaker than the constraint from D^-D^ 
mixing. 

Combining the Mixing and Decay Relations: Assuming that Z' saturates the observed ex- 
perimental value for xd, the bound obtained from the D° — » yU^/x" branching fraction as a 
function of Mz' is 

SGprnlM^XD 



^ Am] 1 



16v^7r5z)r(me, Mz') \ Ml M|, 
2.4x10"^° (sD/M|,(TeV)) < 2.4 x lO^^V^^I'lTeV) . (40) 



C. Family (Horizontal) Symmetries 

The gauge sector in the Standard Model has a large global symmetry which is broken by 



the Higgs interaction 29|]. By enlarging the Higgs sector, some subgroup of this symmetry 



can be imposed on the full SM lagrangian and break the symmetry spontaneously. This 
family symmetry can be global as well as gauged 3^. If the new gauge couplings are 
very weak or the gauge boson masses are large, the difference between a gauged or global 
symmetry is rather difficult to distinguish in practice. In general there would be FCNC 
effects from both the gauge and scalar sectors. Here we consider the gauge contributions. 

Consider the group SU{2)g acting only on the first two left-handed families (it may be 
regarded as a subgroup of an SU{3)g, which is broken). Spontaneous breaking of SU{2)g 
makes the gauge bosons Gi massive. For simplicity we assume that after symmetry breaking 
the gauge boson mass matrix is diagonal to a good approximation in which case Gj^ are 
physical eigenstates and any mixing between them is neglected. Leaving further discussion 
to Refs. [2^, we write down the couplings in the fermion mass basis as 

T<hs = -/ G'i^|sin26'd (rfLT^rf^ - sil^^SLj + sin26'„ {ulIi^ul - cl1^,cl) 
+ sin26'; {cLli^eL - JIlIi^I^l) + cos 26*^ {diliiSi + SLlt^dL^ 
+ cos26'„ {uLlt^CL + CL'jf.UL) + COS 26*, {eilf^lJ^L + T^iT/^ei)} 

+ Gs^jcos {dhltidL - SLltiSLj + cos 26'„ {ul1,,ul - cl^^^cl) 
+ COS 261 {eilf^eL - JIlIi^I^l) - sin 26^ (dilf^SL + SLlf^di) 

- sin 2eu [ulIi^cl + CLlf^ui) - sin 2di {'El'J^.Hl + /^iT/^ei) } . (41) 

12 



Applications of this general interaction yield expressions for D^-D^ mixing, 

(Fs) _ 2flMDBDr{m,,M) ^ ( cos^ ^ sin=^ 2g» 1 \ 

" 3f^ ^ [~^^^^ ' ^^^^ 

for ji^ decay, 



„(FS) _ Mpf^m^ f sin 26'„ cos 261^ cos 26^ sin 201 \ , . 

^DO^u+u- ~ CA^r J Z:2 ZZ2 ' \^^) 



'^Ji r4 I '^^u COS 2U£ COS 2Uu sm 'ZU£ ] ^ 

and for D° /i+e^ decay, 

r(fs) _ Mpflml ( cos29u cos2g^ ^ 1 sin2e^^sin2^y . 

In Eq. P2|) the very-slowly varying RG factor r(mc, M) is set to the scale M ~ 1 TeV. 

Precise predictions for the above three processes are not immediate due to the large 
number of NP parameters. Different patterns can be obtained depending on the region of 
parameter space: 

Case A [mi = -C m2 and 9u — Oe = n/4]: 
Here, B^^o\^^+^- is suppressed and a parameter-free prediction for B^o^^+fj,- in terms of 
occurs. 

Note that here we related the Bdo^^+^i- to the square oi xd- 

Case B [mi = m2 = and 9u — Oe = 7r/2]: 
In this case, the amplitudes for all three processes vanish. 

Case C [mi = m2 <^ m^ and 6^ — = 'n'/2]: 
Now, the mixing contribution vanishes but the branching fractions for fi^fi^ and 

fi^e~ are equal, although undetermined due to NP parameter dependence. 

Case D [mi = m^ ^ 7/12]: 
In this limit, is negligible and there is a parameter-free prediction for Bdo^^+^- 

in terms of x^^, but Xd has the wrong sign. 

IV. NP MODELS WITH LOOP AMPLITUDES 

Although tree amplitudes represent the most obvious situation for producing a correlation 
between mixing and decay, it turns out that loop amplitudes can have the same effect. 
As is well known 3l|, low energy effective lagrangians continue to provide the most useful 



description. In the following, we consider three examples of NP models with loop amplitudes. 
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FIG. 2: Box contribution from heavy weak-isosinglet quarks. 



A. Heavy Vector-like Quarks: Q = -1/3 Singlet Quark 



We first consider models with a heavy vector- hke Q = —1/3 singlet quark. Note that 
essentially identical results hold for a SM fourth quark generation as well, since in each case, 
the fermions will interact with a SM gauge boson and thus the charged leptons have 
SM interactions 32|. It is this which allows for correlations between D^-D^ mixing and 



decay. 

For the class of models with Q = —1/3 down- type singlet quarks, the down quark mass 
matrix is a 4 x 4 array if there is just one heavy singlet (or 6x6 for three heavy singlets 
as in Eq models). As a consequence, the standard 3x3 CKM matrix is no longer unitary. 
Moreover, the weak charged current will now contain terms that couple up-quarks to the 
heavy singlet quarks. For three heavy singlets, we have 



(46) 



where Ui^i = {u, c, t) l and Da = {D, S, B) refer to the standard up quark and heavy isosinglet 
down quark sectors. The {Via} are elements of a 3 x 6 matrix, which is the product of the 
3x3 and 6x6 unitary matrices that diagonalize the Q = +2/3 and Q = —1/3 quark sectors, 
respectively. 

D^-D^ Mixing: The box diagram contribution to AM^ from these new quarks is displayed 
in Fig. [2J Assuming that the contribution of one of the heavy quarks (say the S quark, of 
mass ms) dominates, one can write an expression for x-q j32|. 



(-1/3)1 
D I 



67r2r 



D 



-BD{V:sVusr r{m,,Mw)\E{xs) 



where xs 



[ms/MwY . The Inami-Lin 



E{xs) = Xs 



1 



9 



3l[ | function E{xs) is defined as 
3 



+ 



3x| 



4 A{xs-l) 2(x5-l)2 2{xs-l) 



Inxs 



(47) 



(48) 
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For our numerical work, we assume a default value of ms = 500 GeV, but express our 
result for variable ms by noting that the functions E{xs) and C{xs) (c/Eq. fl53|) below) are 
proportional to xs within ten per cent over the mass region 400 < ^^^(GeV) < 700. The 
light-heavy mixing angles iK^V^usI^ should go as l/ms for large ms to keep the contribution 



under control. The current bound on 



very stringent, \V*sVusf < 4 x lO"'' [l6|. 



2 

^5^5 1 froKi unitarity of the CKM matrix is not 



In the -Eg-based model proposed by Bjorken et al [33|], the 6x6 mass matrix has an 
especially simple form. The resulting 6x6 mass matrix has a pseudo-orthogonality property 
which implies that the 3x3 CKM matrix, although not unitary, satisfies 



(49) 



i=X 



The analog of this condition in the up quark sector does not hold, and as a result, there 
are no new FCNC effects in the down quark sector. For the CKM elements participating 
in D^-D^ mixing, the prediction is now (recall capital lettering is used to denote the heavy 
quark) 

(50) 



KsVuS^" = sl\V:,Vus\^ sW 



where IV^^^msI ~ A ^ 0.22 and S2 is the (small) mixing parameter describing the mixing 
between the light s quark and the heavy S quark. Thus, we rewrite |2;d"^''^''| in the modified 
form, 

(-1/3) I ^, G^M^flMo ^ ^2\2, 



\x 



67r2r 



D 



-BDsiyr{m,,M^)\E{xs)\ . 



(51) 



D° — ^ 11^ II Decay: For fi^fi , the effective Lagrangian is given in Eq. (2.2) of 



Ref. 



3l| 



eff 



where 



C{xs) ^ 



7r2 



3x,c 



C{xs)\S2 Ql , 



(52) 



Xq 



4 4(x5 - 1) 4 \xs - 1 
In this model, the — branching fraction becomes 



\nxs 



(53) 



Mr 



Bdo- 



nJl - 4^ {GfM^)^ ■ {s2\fDm^C{xs) 



(54) 
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Combining the Mixing and Decay Relations: If we eliminate from the mixing and decay 
relations, we obtain 



6 xnJl-^ml/Ml (m^GFMwC{xs) 

Br)0 — 



-l.OxlO-^xo ( y < 1.0 X 10-- ( y . (55) 

V500 GeVy - V500 GeVj ^ ' 

B. Minimal Supersymmetric Standard Model 

We next consider the Minimal Supersymmetric Standard Model (MSSM) with unbroken 
R-parity. Conservation of R-parity implies that only pairs of sparticles can be produced or 
exchanged in loops. We will assume that neither squarks nor gluinos are decoupled (direct 
collider searches for squark and gluino pair production place the bound m^ ^ > 330 GeV [l6| 
in the MSSM with minimal gravity mediated Supersymmetry breaking), so the MSSM can 
in principle give a dominant contribution to the processes under consideration here. 

We will not assume any particular SUSY breaking mechanism, and hence parameterize 
all possible soft SUSY-breaking terms. We work in the so-called super-CKM basis, where 
flavor violation is driven by non-diagonal squark mass insertions (see [l^ for a discussion 
of this mechanism in — mixing). In this case, the squark-quark-gluino couplings are 
flavor conserving, while the squark propagators are expanded to include the non-diagonal 
mass terms. The 6x6 mass matrix for the Q = +2/3 squarks can be divided into 3x3 
sub-matrices, 



and the mass insertions can be parameterized in a model independent fashion as 



(56) 



Here, i,j are flavor indices, M,N refers to the helicity choices LL, LR, RR, and mq repre- 
sents the average squark mass. The squark-gluino loops with mass insertions are by far the 
largest supersymmetric contribution to D° — i5°-mixing and can dominate the transition. 
The effective hamiltonian relevant for this contribution to — D°-mixing is given by 

^MssM = 7r^j:Ciim,)Q, , (58) 

^"'g i=l 
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where all eight operators contribute in the MSSM. Evaluating the Wilson coefficients at the 
SUSY scale gives, 







= ^(^r2)L [4a;/i(x) + ll/2( 








C2 




= ^{W2UW2U 15/2(0;: 


-(5r2)LL(^r2)«« 


[2x/i(x) + 10/2 






H) 


= lm2)LL (5r2)«K[42x/i(x 


) - 6/2(x)] - (5^2) 


Lij('^r2)iJL ll/2( 


x)} 




H) 


= l^(^r2)k37xA(x), 








c, 




= ^(^r2)k^/i(^)' 










H) 


= ^(^r2)L [4x/i(x) + ll/2 








C7 




= ^(^r2)L37x/i(x), 








Cs 




= ^(^r2)L ^/i(^), 









(59) 



where x = m'^/m'^, with being the mass of the gluino. The equations above a re sy mmetric 



under the interchange L ^ R. These contributions to x^i are found to be large [l^l, and the 
observation of D mixing constrains the mass insertions to be at the percent level, or less, 
for Tev-scale sparticles. 

We now examine the squark-gluino contribution to rare decays, which proceeds through 
Z penguin diagrams for on-shell leptons. The relevant c — ul^l~ Lagrangian is given, 
for example, in Ref. 20|. Electromagnetic current conservation forbids the contribution 
of the photonic penguin diagram for on-shell leptons in the final state. In addition, the 
vector leptonic operator C'-y^l also does not contribute for on-shell leptons as p^i^'^fii) = 
{p^+ +p^_)(£7^£) = 0. The effective hamiltonian is then given by 

A 2 

where Cio and c[q are given by jsJ ] 



1 rv 1 
ClO = -^-{5^22)LR{5'i2)RLP032 , c[, = ^^(5J^^)^^{S-^^)^^P,,, . (61) 

9 a 9 a 

-Po32,i22 are kinematic loop functions and are defined in the above reference. The double 
mass insertion is required to induce a helicity flip in the squark propagator. Due to this 
double mass insertion, this contribution to D ^ is completely negligible. We note 

that the chargino contribution to the Z penguin for D £^i~ also contains a double mass 
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'Li 



c u 

u c 



FIG. 3: Contributions to D^-D^ mixing from the A' superpotential terms in super symmetric models 
with R-parity violation. 

insertion. The leading MSSM contribution to this rare decay is thus most likely mediated 
by a box diagram with squark-chargino-sneutrino exchange. This precludes a relation to 
— mixing. 

Lastly, we note that in contrast to the Bs system [35|, D — > does not receive a 

sizable contribution from Higgs boson exchange with large tan/5. This is because in this 
case, the loop-induced term to the Yukawa couplings is proportional to Vd {i.e., the vev of 
the Higgs doublet that generates masses for the down-type quarks) which is the smaller of 
the two vevs and hence does not compensate for the small loop factor. 



C. R Parity Violating Supersymmetry 



Finally, we consider Supersymmetry with R-Parity violation (RPV). We refer the reader 



to Refs. 
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20| for discussions and earlier references of RPV-SUSY relevant to this paper. 



Suffice it to say that the lepton number violating RPV-SUSY interactions can be expressed 



as 



(62) 

in terms of the coupling parameters {A^^j^,}. The generation indices denote the correspon- 
dences i -x^ leptons or sleptons, j up-type quarks and k down-type quarks or squarks. 

D^-D^ Mixing: As described at the high mass scale by the effective hamiltonian 



12871^ 



1 1 



2 



Qi , 



(63) 



this implies constraints on the product of couplings X'i2kKik- Here, we have assumed that only 
one set of the R-parity violating couplings X'i2kKik large and dominant. This is equivalent 
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to saying that, e.g., both the sleptons and both the down-type quarks being exchanged in the 
first contribution to the box diagram of Fig. [3] are from the same generation. In general, this 
need not be the case and the couphng factor would then be the product K2k^'mik^'m2nKiny 
with, e.g., the set of iL,i,dR,k,(^L,m,dR,n being exchanged. Computing the evolution to the 
charm-quark scale yields 



n 



Rr, 



2m^ 



-Ci{mc)Qi 



(64) 



with Ci(mc) = r{mc,mq)Ci{mq). The mixing contribution from the R-parity violating A' 
terms then implies 

1927r2rDm2- 



where e = m'] /m? 
300 GeV, so that 



iK2kKlk)' 



{l + e)flMDBDr{mc,mq) 
For definiteness, we shall scale the results to the value mi 



(65) 



0.0053 ■ 



300 GeV 



1 + e 



Xb 
0.01 



(66) 



fj^~^fi~ Decay: In RPV-SUSY, the underlying transition for fi~^fi~ is c + m 

fi^ + /i~ via tree-level d-squark exchange. The coupling constant dependence for D'^ - 
would therefore generally involve A^a/c-^iifc- For the specific mode D° /^^A*^, we take i = 2 

i ^ fi, reads 

(67) 



to get X'22k^2ik- The effective hamiltonian of Eq. (74) in Ref. (20|], but with 



ortes — — 9 — VI 



This leads to the branching fraction 



B 



DO- 



D 



2^1 



4m?, 



Ml 



1/2 



('^22A:'^21A: 

1287r m^- 



and so the constraint 



(-^22fc-^21fc) < ^D^'- 



1287rmS r 

dk 



D 



foK Mb 



Mr 



2 1-1/2 



which reads numerically 



(68) 



(69) 



AmA'sifc < 0.088 



("300 GeVy 



dk 



(70) 
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Model 




Experiment 


< 1.3 X 10^^ 


Standard Model (SD) 


~ 10-18 


Standard Model (LD) 


^ several x IQ--*^^ 


Q = +2/3 Vectorlike Singlet 


4.3 X 10-11 


Q = —1/3 Vector like Singlet 


1 X 10-11 (m^/SOO GeV)2 


Q = —1/3 Fourth Family 


1 X 10-11 (mcj/SOO GeV)2 


Z' Standard Model (LD) 


2.4 X 10-12/(^^2' (TeV))2 


Family Symmetry 


0.7 X 10-18 (Case A) 


RPV-SUSY 


4.8 X 10-9 (300 GeV/mjj2 



TABLE L Predictions for n^fi branching fraction for ~ 1%- Experimental upper bound 

is a compilation from 



Combining the Mixing and Decay Relations: Note the mixing constraint involves A^gfc-^iifc) 
whereas the decay constraint has y22k'^2ik- i = 2 case dominates, then we arrive at 

the prediction (here we set e = 1) 



'^DO^fj.+fi- 



1 - Aml/MY) 



1/2 



Xd 



4m| BDr{mc,mg) 



4.8 X 10~^ Xd 



'300 GeV' 



< 4.8 X 10" 



'300 GeV' 



mi 



m 



(71) 



V. CONCLUSIONS 



The search for New Physics will in general involve many experiments, including the mea- 
surement of rare decay branching fractions and observation of particle-antiparticle mixing. 
Such experiments are essentially competitors, each seeking to be the first to indirectly de- 
tect physics beyond the Standard Model. At any given point, which measurements are more 
sensitive to New Physics must be determined on a case by case basis. Our earlier work of 
Ref. already pointed out that the observed D^-D^ signal imposes severe limits for a 

20 



large number of New Physics models. If the D^-D^ mixing is dominated by one of those 
New Physics contributions, what does it imply for rare decays such as —>■ /i+yU"? Not 
only have we been able to answer this question in several specific scenarios, but we find 
a striking correlation in some of the models, wherein the branching fraction for the decay 
mode fi'^fi~ is completely fixed in terms of the mixing parameter Xd- 

For convenience we have gathered our results in Table [B For all but one case (Family 
Symmetry), we find the NP branching fraction exceeds the SM branching fraction. All the 
NP branching fractions are, however, well below the current experimental bounds of Eq. ([2]). 
Anticipating future improvements in sensitivity, the first NP model to be constrained will 
be R-parity Violating supersymmetry. This will require lowering of the current bound by a 
factor of a few hundred. 
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APPENDIX A: RG RUNNING AND MIXING MATRIX ELEMENTS 



NP contributions are affected by RG-running. The relevant 8x8 anomalous dimension 



matrix, derived to NLO in Ref. 
Wilson coefficients, 



2l| . is applied to LO here (see also Ref. [l^) to yield the 



Ci(/i)=r(/i,M)Ci(M) 
C2(/i)=r(/i,M)i/2^2(M) 
C3(/i)=r(^,M)V^^^^ + r(/.,M)- 



C4(/i)=r(^,M)(i+^)/^ 



V24T 



Cs{M) 



Ci{M) 



2C2(M) 



30C5(M) 
V24I 
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C5(/i)=r(/i,M)(i+v^)/^ 
where (presuming that M > rrit), 



i 8 
2^ 



V24T 



C4(M) + 
C5(M) + 



30C5(M) 

Ci{My 

C4(M) 
8^241 



(Al) 



r(/i, M) 



(A2) 



/ Us{mt) \ ^^^^ / as(m5) \ ^^^^ 
as{mt)) \as{mb) J \ a,{fi) J 

Regarding the remaining Wilson coefficients, Cq runs analogous to Ci and Cy § run analogous 
to C4,5. The presence of operator mixing in Eq. (lAip is a consequence of the nondiagonal 
structure of the anomalous dimension matrix. 

We also need to evaluate the D^-to-D^ matrix elements of the eight dimension-six basis 
operators. In general, this implies eight non-perturbative parameters that would have to 
be evaluated by means of QCD sum rules or on the lattice. We choose those parameters 
(denoted by {Bi}) as follows. 



{Qi) = 


IflMlB, , 




= flMlB, , 


{Q2) = 


-yiMiB, , 


m 


= IflMlB, , 


{Q3) = 




{Q7) 


= -yiMiB 


{Qa) = 




{Qs) 





(A3) 



where {Qi) = (-D°|Q,;|-D°), and represents the D meson decay constant. By and large, the 
compensatory 5-factors {Bi} are unknown, except in vacuum saturation and in the heavy 
quark limit; there, one has Bi —>■ 1. 

Since most of the matrix elements in Eq. (1A3|) are not known, we will need something 
more manageable in order to obtain numerical results. The usual approach to computing 
matrix elements is to employ the vacuum saturation approximation. However, because 
some of the 5-parameters are known, we introduce a 'modified vacuum saturation' (MVS), 
where all matrix elements in Eq. (1A3P are written in terms of (known) matrix elements of 
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{V - A) X {V - A) and {S - P) x {S + P) matrix elements 5d and b'^ 



(S) 



(Qi) = ^flMlBn 
{Q2) = flMlBo 



1 

"2 



1 7] 

4iVe^ 2 
f!,M^Bn V 



(Q5) = yJImIBo V 

m = {Qi) , 
iQi) = m , 

{Q8) = {Q5) , 



where we take iVc = 3 as the number of colors and define 



V 



< MB 
Bn m?.. 



In our numerical work, we take 



(A4) 



(A5) 



1. i?D = 0.82, which is the most recent result from the quenched lattice calculation. 

2. For 7], we use B^^ ~ Bd [36] so that r] ~ M'^/ml ~ 2. 



3. Regarding the decay constant there is now good agreement [37[ between determi- 
nations from QCD-lattice simulations /^^**''' = 0.207(4) GeV and various experiments 
^^expt.) ^ 0.206(9) GeV. For definiteness, we adopt the value fo = 0.207 GeV. 
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